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1 Introduction 



Dynamical systems come with a variety of behaviours, which have been the object of 
many studies. In particular, the long term predictability is a subject of interest since 
the first considerations on the solar system and received their modern formulation from 
Poincare a century ago 

It has been shown that birational maps have a great variety of dynamical be- 
haviours, in spite of their apparent simplicity: We need only the basic four operations 
to compute the successive points and we can even dispense with the division if we 
keep with homogeneous coordinates. The purpose of this note is to relate different 
characterizations of dynamical systems for birational maps. 

The first one is the existence of invariants, which allow for a Liouville like form 
of integrability. An other one is the sensitivity to initial conditions, which is related 
to the derivative. The last one is the "algebraic entropy" [^. Common factors in 
the homogeneous coordinates of the image of a point result in simplifications of its 
expression and decrease the degree of the n-th iterate of the transformation from the 

value obtained without simplifications. This notion is directly linked with Arnold's 
complexity [Q], since the degree of a map gives the number of intersections of the 
image of a line and a hyperplane. 

A zero value of the algebraic entropy corresponds to a sub-exponential growth of 
the successive degrees and in all known cases to a degree depending polynomially on 
n. A first proof of such a polynomial bound on the degrees was given in [^]. It has 
been observed on numerous examples through the determination of the images of a 
generic line. This would result from the conjectured existence of recurrence relations 
on the degrees. 

This work is motivated by the observation that whenever the iterates are confined to 
a low dimensional subvarieties by conserved quantities, it has been possible to find that 
the algebraic entropy is zero. Moreover, in many cases the growth is simply quadratic. 

Here I show how this polynomial behaviour of the degrees can be derived from the 
existence of invariant subvarieties which are abelian varieties. In the simplest case of 
rational or elliptic curves, the exact exponent will be proven. 

The first section will precise what is called integrability in this context of discrete 
time dynamics. Next will come an analytic derivation of a quite general bound on the 
growth of the degrees of the iterates, based on the expression of the degree of a variety 
from its volume. The last section makes use of the notion of addition on an elliptic 
curve to build iterates with degrees of quadratic polynomial growth. The practical use 
of such constructions and the foreseeable extensions are the subject of the conclusion. 

2 Action-angle variables 

Even if it is intuitively clear that integrability has to do with the existence of a sufficient 
number of conserved quantities, the question remains of what is sufficient. A discrete 
time analogue of Hamiltonian mechanics has been proposed [^], allowing for a direct 
generalizations of the answer provided by Liouville [^. However, the resulting class 
of systems is too restrictive for our porpoise and we cannot rely on a definitive answer. 
Since even two-dimensional algebraic varieties can have birational maps with com- 
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plex dynamical behaviour, the only safe case is the one of one-dimensional invariant 
varieties. In this case, the induced birational map on the invariant varieties are in fact 
holomorphic: Singular varieties of birational maps are of codimension at least two 
and there are no such subvarieties in curves. It remains to classify automorphisms of 
infinite order of the curves, which exist only for curves of genus or 1. Higher genus 
curves have but a finite number of automorphisms. The only possible automorphisms 
have therefore finite order and as dynamical systems, they are periodic. 

In the genus case, the curve is the Riemann sphere and the automorphisms are 
homographic transformations. By a reparameterization, we can always obtain one of 
the cases g{z) = z + a or g{z) = \z, according to the number of fixed points of the 
transformation. In the multiplicative case, we can describe it as a translation if we set 
z = exp t. The point depends on the parameter t with the period 27r/ and now the 
transformation is a simple addition on t. 

Non-singular genus 1 curves are holomorphically isomorphic to a torus C/Z + rZ. 
There is therefore a parameterization of these curves by doubly periodic functions of 
z G C. In terms of this parameter, the infinite order automorphisms are translations 
z ^ z + 13. 

In all these cases, the map is described by a simple translation in a parameter 
which can be one to one, simply periodic or doubly periodic. The n-th iterate of the 
map can be simply expressed as the translation by the same parameter multiplied by n. 
However there is a price to pay, especially in the genus 1 case, since we have no longer 
algebraic functions, but a parameterization by WeierstraB functions or other elliptic 
functions. 

Nevertheless, the scheme is that in adapted coordinates, the transformation takes 
the simple form: 

ct>{x,e) = {x,e + a{x)), (1) 

with X designing the invariants of the integral curve and 6 a parameterization of this 
curve. In these coordinates, the n-th iterates of is easy to write: 

<P''{X,e) = {X,e + na{X)). (2) 

From the explicit form of (/>", its differential is deduced to be: 

The important property is that the matrix elements grow at most linearly in n, giv- 
ing a computable behaviour for long times, since numerical errors get multiplied by 
small numbers. In generic systems, there would be hyperbolic fixed points of cj)^ . At 
these points, the differential of ij)^^ grows exponentially. Small errors in the first few 
iterations can have sizable effects after a moderate number of iterations. 

There are natural generalizations of this scheme for higher dimensional varieties, 
but in the absence of an equivalent to the symplectic structure which is essential in the 
theorem of Liouville, they will not be necessary outcomes of the existence of invari- 
ant varieties. Even in a two-dimensional variety with rational parameterization, the 
restriction of a birational map can have complex behaviours. And it is not even suffi- 
cient that the map restricted to the invariant varieties has no singularities, since a torus 
of two complex dimensions can have exponentially diverging trajectories. An example 
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is the following bijective transformation of /I? + tI? (the product of two elliptic 
surfaces with the same modulus): 



(01,62) ^{291 + 02,61 + 92). (4) 

The following section is equally valid for the case of abelian invariant varieties if 
the restriction of the map to these varieties is a simple translation. The formulas 
apply with a multidimensional 6. Examples of such systems were shown to exist in [Q], 
with a map which reduces to an addition on the Jacobian of a curve. It is however not 
clear if these systems can be expressed as maps on projective spaces, since they are a 
priori defined on the quotient of a linear space of matrices by gauge transformations. 
In the case which gives one-dimensional Jacobians, it is possible to fix the gauge and 
the system reduces to the one associated to the height- vertex model of Baxter [^. In 
the other cases, it is not clear whether such a choice is possible in general, without 
restricting to a lower dimensional part of the Jacobian. Part of the symmetry group of 
a three dimensional model of statistical mechanics gives a dynamical system in CP^ 



with two-dimensional invariant varieties []10|]. From the images of the orbits, it seems 



that the restriction of the map on an invariant variety is a shift on a torus of complex 
dimension two, but this has yet to be proved. 



3 Analytic method. 
3.1 Fixed points 

An immediate consequence of the form of the compositions of (/> is that the fixed 
points of the iterates are not isolated, but form whole varieties. If p is a fixed point of 
(/)", it means that na{X) is one of the periods of the pai-ameterization of the invariant 
curve it belongs to and all points of this curve are fixed points of cj/^. Moreover, the 
equation determining the fixed point variety reduces to an{X) = 0. This is as many 
equations as there are components in 6. The codimension of the fixed point varieties 
is therefore the dimension of the invariant varieties. 

In the bidimensional case, the fixed point variety is of the same dimension as the 
invariant varieties and reduces to a product of invariant varieties. This factorization is 
however not straightforward since the corresponding values of the invariants are gener- 
ically irrational algebraic numbers. For the determination of the individual invariant 
varieties, the factorization of the equation of the fixed point variety must be done in an 
algebraic extension of Q which has to be determined. 

In higher dimensions, the fixed point variety would correspond to a variety in the 
space of invariants, so that it would be still more difficult to deduce the equations of 
individual invariant varieties. However, the main problem in determining invariant 
varieties is to find their covariance factor, that is the factor they get multiplied by when 
applying the transformation. The fixed point varieties are invariant and their covariance 
factors are powers of the covariance factors of the invariant varieties. 

Moreover, the study of the differential of the map at fixed points should be a good 
starting point for a proof of the existence of invariant varieties. 
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3.2 Analytic proof of a polynomial growth of the degrees 



The degree d of an algebraic subvariety is proportional to its volume in the projec- 
tively invariant Kahler metric of CP^ [11]. The basic idea is that the volume forms 



induced on complex subvarieties by a Kahler metric are simply exterior powers of the 
associated Kahler form. The Kahler form being closed, this volume is invariant by 
deformations and in particular, can be computed from the retraction of the variety to 
a linear space L of the same dimension. In a generic situation, points outside a lower 
dimensional subspace of L have d preimages by the retraction, so the volume of a 
variety is d time the volume of a linear space. Since all linear spaces are related by 
U {N + 1) transformations which leave the metric invariant, the volume of a linear 
space is a universal constant. 

This gives an analytic way of computing the degrees. The degree of the map (j) is 
the degree of the image of a generic complex line L and this degree can be computed 
by an integration, since it is proportional to the area S: 

S = [ oj= [ r.u:. (5) 

J(/>"(L) JL 

In this formula u> denotes the Kahler form. The definition of its pull-back (p^uj involves 
the differential dcj)"- of 0" to transform the tangent vectors of L into tangent vectors 
of its image. More precisely, the form uj and the differential dcj) are sections of fiber 
bundles with values at a point x denoted respectively by uj^ and d4>x- is a bilinear 
form on the tangent space T^M, d(pa is a linear map from T^M to T^(^i£^M. (p^u) is 
given by 

{ct),u)a{X,Y) = UJ^^,){d^aX,dcl)aY), (6) 

with X, Y in T^M. From the formula (js]) for d^" follows that we have a simple de- 
pendence on n of the integral in (js]). The Kahler form must be evaluated in the variable 
point 0"(a), with the tangent vector expressed on a basis associated to the action-angle 
coordinates. CP^ is compact, uj is bounded in any affine coordinate patch from its 
explicit expression and the change from action-angle variables to projective coordi- 
nates is smooth, so that to remains bounded when expressed on an adapted basis of the 
tangent space. The boundedness of lOx is then sufficient to conclude that the degree of 
(i" is less than some constant times n?. 



4 Algebraic proof 
4.1 Generalities 

In the preceding section, I showed limits on the degrees of the iterated maps. The proof 
however does not show how to build transformations of the given degree. This will be 
remedied in this section. 

The basic idea is that eq. can be expressed algebraicly, without the introduction 
of a specific parameterization of the invariant variety. Since every analytic operation 
on an algebraic variety is algebraic, A + B is an algebraic function of A, B and 
parameters describing the invariant variety. The algebraic translation of eq. @ is made 
of operations independent of n and the computation of na. The whole n dependence 
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comes from this operation and we want to find a good bound of the degree of this 
operation. All other operations can contribute an overall factor or some additional 
constants but cannot modify the growth behaviour. Computing 2k X from /c X is a 
simple addition. 2^ X can therefore be computed by p additions and for a general n, 
the number of additions to perform to calculate n X is of order log2 n. If the degree 
of the result is multiplied by a constant r by each addition, the degree of n cr will be 
of order r^°S2 which is equal to n^°^'^ ^ . We therefore have built an expression for the 
image of a point with a degree satisfying the required bound. 

In fact, this discussion is somehow naive. The invariant variety depends on the 
starting point and care must be taken of the dependence of the addition on the param- 
eters of the invariant variety, which will in turn be some polynomial expressions of 
the coordinates of the starting point. Going from the point pz to the point 2pz is an 
operation of degree d in the coordinates of the point pz and of some degree I in the 
coordinates of the initial point. Let us call nik the degree of the operation 
nij^ satisfies the following relation: 

mfc+i = druk + I. (7) 

It is elementary to verify that is also bounded by Cd^, except in the trivial case 
d=l. 

The degree of the iterate of order n of an integrable mapping is therefore bounded 
by a polynomial function of n. In the following, the precise form of this bound will be 
established in the case of curves. A generalization for higher dimensional cases is in 
principle straightforward, but rather involved. 

4.2 Curves 

The case of rational invariant curves can be easily settled, since the parameterizations 
of the invariant curves are rational. A birational change of variable allows for a product 
structure CPi times the space of the invariants. 

In the case where the ^-variable is directly the rational parameter, the expression 
of 0" can be directly read from @. In this case, the successive iterates have all the 
same degree. 

In the other case, the transformation takes the form t — > \{X)t. The iterates 
depend on \{X)^ and the degree of the transformation is linear in n. 

In the case of elliptic invariant curves, the choice of a special point O allows for 
the identification of the curve with its Jacobian. This determines the addition on the 
curve: the sum S of the points A and B is the second zero of a meromorphic function 
with simple poles at these two points and a zero at the origin O. Such a meromorphic 
function can be obtained as the ratio of sections of a line bundle. In the simplest cases, 
these sections are simply linear functions of the homogeneous coordinates. Their ze- 
roes are given by the intersection of the curve and some hyperplane. 

There are two simple algebraic descriptions of an elliptic curve. The first is a de- 
gree 3 plane curve, the second is the intersection of two quadrics in three-dimensional 
space. In both cases, the explicit determination of multiples of a point will allow for a 
quadratic bound on the degrees of the successive iterates. 

Otherwise, we could parameterize the curves by doubly periodic functions, Weier- 
straB function for the degree 3 plane curve or ^-functions for the biquadratic. Then 
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addition formulas for these functions would translate in algebraic formulas for the de- 
termination of the sums, but I do not want to introduce any considerations on these 
special functions. 

Higher degree representations are singular with multiple points. These multiple 
points correspond to multiple values of the parameter of the curve and are sent to 
a number of differing points. They are therefore singular points of the transforma- 
tion and it is possible to make birational transformations which resolve those multiple 
points without introducing other multiple points. 

4.3 Degree 3 plane curve 

Let us consider a curve with the equation: 

V^ = Pi{x). (8) 

P3 is a degree three polynomial with no multiple roots. The natural choice for the 
origin O in this case is the point at infinity. This point is common to all invariant 
curves so that the parameter cr{X) cannot be obtained directly as the image of the 
origin O, but must be computed as 0(P) — P for a generic point P on the curve. 

Affine functions of x and y are sections of a line bundle when restricted to the 
curve. The zeros of these sections are given by the intersection points of a line and the 
elliptic curve: there are three of them. 

To build a meromorphic function with poles in A and B, we first consider the line 
through A and B: it will cross the curve in a third point M. We will take now a line 
going through M and the origin O: I claim that the third intersection point of this line 
with the curve is the sum S. Indeed, taking the ratio of the two above mentioned linear 
functions, we get a meromorphic function on the curve with poles in A and B and 
zeros in O and S. The zeros in M of the two linear functions cancel each other. Lines 
going through O have equation x = xq and intersect the curve in two points of the 
form (x, ±y). 

When defining the double of a point, we need to parameterize the tangent to the 
curve at the given point. It will be given by (x + 2ay, y + aP^{x)). Assuming a 
reduced form for P3, that is P3 = x'^ + ax + b, the equation for a to get a point on the 
curve reduces to: 

P^(x)2a2 = Sa^y^ + Ua^xy"^. (9) 

The non zero solution for a can be substituted and expressed in terms of x from the 
equation Remembering that the point 2z has the opposite y from the intersection, 
we obtain the following coordinates for it: 

/ P^(x)2 / 3xP^(x) Piixf \\ 

The important point is here that the new x is of degree 4 in the old x, without any 
dependence on y. As in the preceding section, it is possible to deduce that the degree 
of X for the A;-th iterates is at most quadratic in k. The degree of y does not seem to 
be as easily bounded. But the very expression of the new y shows that its degree is the 
degree of the old y plus 6 times the degree of the old x. The degree of the variable 
y after a number of steps will be 1 plus some times the sum of the degrees of the 
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preceding x. This is the sum of a geometric series, so that it will also be of order 4^, 
apart from some constant prefactor. 

In this case, the integer d is 4, that is to say the bound derived in the preceding 
paragraph is a quadratic one. 

4.4 Biquadratics in 3 dimensional space. 

The other possible equation for a non-singular elUptic curve is that of a biquadratic in 
a three dimensional space. By a proper choice of coordinates and of combination of 
the equations, they can be brought to the form: 

— ay^ — hz^ = 0, 
t^-hy^-az^ = 0. (11) 

This readily shows the symmetries of the curve, since these two equations are invariant 
by the change of sign of any of the coordinates. Since the change of the sign of 
all coordinates is the identity in projective space, this gives a group Z2 x Z2 x Z2 
of symmetries. Section of the tautological bundle are just given by a linear form in 
x,y,z,t. They are defined up to a factor by a plane in this three-dimensional space, 
that is by three points. Intersection of a plane with the curve will generically consist 
in 4 points. Computation in the Jacobian can easily be done if we take as zero a point 
with one of the coordinates 0. There are four of these points on any curve, but they are 
related by the afore mentioned symmetries. 

A plane is defined by three points, so that we choose a point R on the curve which 
will belong to the planes corresponding to the two factors defining a meromorphic 
function. The plane defined by A, B and R will have a fourth intersection with the 
curve, H, and the plane going through H, R and O will meet the curve at the sum 
S: when taking the quotient, the zeroes in H and R cancel and we get the required 
meromorphic function with the divisor A + B — S — O. 

R should be chosen to be one of the points with the same zero as the origin, since 
in this case S is simply related to H. In this case, the direction of the Une OR is one 
of the coordinate axis. It is either the axis where the coordinates of O and R differ 
by a sign or the axis of the null coordinate in the case R and O coincide, since the 
tangent in this point is simply giving by the vector whose only non-zero component is 
on the null coordinate. Then if we take for S the point obtained from H by negating 
the coordinate along the axis OR, HS and OR have the same direction and they are 
coplanar as two parallel lines define a plane. A remark which be useful in the sequel 
is that if O is the point (xq, yoj zq, 0) adding the point Oz = (xq, yo) — -^Oj 0) to a point 
is simply changing the sign of the z and t coordinates. With obvious notation, the 
addition of the Ox and Oy points are similarly obtained as symmetries of the curve 
changing the sign of two of the coordinates. Since the negative for the addition on the 
curves is obtained by changing the sign of the t coordinate, all obvious symmetries of 
the curve can be simply obtained by a combination of the addition of one of the points 
Oi and the opposite. 

The obvious way of adding multiple times the same point is to define a plane 
which is tangent to the curve at this point, but in this case, we do not end with a result 
of sufficiently low degree. It is however possible to consider the plane going through 
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the original point and two of its transforms by symmetries to compute a triple. The 
plane going through A, A + Ox and A + Oy intersects the curve in the point —3^4 + Oz 

{Ox +0y = OJ. 

If we start from three points Ai, A2 and A^ on the curve, the plane they define is 
parameterized by the linear combination aiAi + a2A2 + a^A-^. When substituting in 
the equations of the curve, the af terms disappear since the Ai 's satisfy these equations 
and we have two linear equations for the three products aia2, aia^ and 0203, whose 
coefficients are bilinear in the coordinates of the points. These products are then of 
degree 4 in the points, the a-i of degree 8 and the final point is of degree 9. We obtain 
a variant of the result we announced. Here, we proceed by tripling and the operation 
is of degree 9, which is the cube of 3, so that we still obtain a quadratic bound for the 
degrees of the iterations of the maps. 

There still remains the problem that in fact, the origin = (xq, yo) zq, 0) of the 
curve is not algebraicly determined. Xq, and Zg are expressed as quadratic polyno- 
mials of the initial point, but the actual coordinates are square roots. We will have to 
keep track of the usage of each of these quantities to be sure that the end result only 
depends on even powers of these coordinates, so that the result is truly an algebraic 
one. Here again, the symmetries of the curve give us the possibility to simply derive 
what happens if we change the sign of one of these coordinates. Changing the sign of 
zo for example, is changing O to O^. Multiplying the image of the origin will give the 
iterated image of this point. Then to obtain the image of the starting point P, we have 
to add P with the image of the origin. But this addition depends on the origin, since 
A + B depends on the fact that A + B — S — O is a. principal divisor. The change 
in the origin completely compensate for the change of the point. The conclusion is 
that the calculation does not depend on the possible choices of signs. Except from a 
global sign factor, the calculated point is then invariant, so that it only depends on even 
powers of the coordinates of the origin and is therefore an algebraic function of the 
starting point. This is an elementary case of Galois theory, which says that element of 
an algebraic extension which are invariant by all elements of the automorphism group 
(the Galois group) are elements of the base field. 



5 Conclusion 

Integrable mappings have been proven to have zero algebraic entropy. Moreover, the 
observed quadratic polynomial growth of the degree is seen to be generic from the 



results of sec. (3.2). 

This step forward in the study of the dynamical behaviour of birational maps urges 
to find a reciprocal result. However the case of the discrete Painleve system studied 
in [^] shows that a zero algebraic entropy does not imply the existence of invariants, 
even in a two-dimensional setting. Maybe this is due to the non-autonomous charac- 
ter of this map, which implies that fixed points of cp'' are not fixed points of all the 
transformations cf)^^. 

The derivation of the higher order iterates sketched in sec. (^ has been useful for an 
alternate proof of the polynomial growth of the degrees. We can wonder whether this 
calculation can be put to more practical uses. The logarithmic growth of the number 
of operations should in any case become a clear advantage both for the precision of the 
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calculation and the required time. 

The usefulness of the fixed point varieties for a characterization of the invariant 
varieties should be probed. But for higher dimensional systems, the explicit equations 
are very large and difficult to obtain with current algebraic manipulation programs. 
It should be of interest to obtain such varieties from implicit equations, i.e., without 
trying to have a formula for the Otherwise, intersection of some low dimensional 
variety with the fixed point varieties should be sufficient to determine their dimension- 
aUty and therefore the dimension of the invariant varieties. 

It would also be interesting to have examples of the unusual behaviours which 
seem possible if the invariant varieties are higher dimensional abelian varieties. As is 
pointed out in eq. (Q), some dynamical exponents could be positive with a zero alge- 
braic entropy, giving "simple" dynamics with diverging trajectories. If the differential 
has a more complex structure than in (|3|), it should also be possible to have a poly- 
nomial growth of the degrees, but faster than quadratic. This would however imply a 
rather long recurrence relation for the degrees, which itself is the sign of a complex 
scheme for the resolution of the singularities 



A Birational maps. 

In this appendix, I want to recall some properties of birational maps which have been 
described in more detail in [Q] and which are relevant to this work. 

Birational maps are not really maps. Except in the simplest case, they are not 



defined everywhere [11]. They are more general correspondences, which can be one 
to many at some points. But they are nevertheless maps to a good approximation since 
they define maps on the complement of an algebraic variety. 

The projective space CP" is the space of complex lines through the origin in C"+^. 
Any non-zero element of the line is a set of homogeneous coordinates for this point 
of projective space. If one of the coordinates is set to 1, the others give the so-called 
affine coordinates of a part of the projective space. 

A rational map between the projective spaces CP™ and CP" is simply a homo- 
geneous polynomial map from C™+^ to C"+^. The homogeneity degree is called the 
degree of the map. The image of the line representing a point in CP™ is either iden- 
tically zero in C"^^ or a complex line by the homogeneity of the map. In the first 
case, the projective point is singular for the map and in the second, the line defines the 
image in CP". Going to affine variables, the map is defined by rational functions, but 
this leads to spurious singularities, which are just the consequence of saying that some 
hyperplane of projective space is at infinity. 

The composition of rational maps requires some care. If we compose the poly- 
nomial maps, (/) of degree and ip of degree give a map of degree d^d^, ip o (f>. 
But with this composition law, only maps of degree 1 can have an inverse, since I is 
of degree 1. However, the greatest common divisor of the coordinates of the image 
can be factored out and the map is not modified at the points where it was defined. It 
is naturally extended to points where this divisor vanishes giving a reduced product 

ip o (j) = m{ip, 4>)iIj X (j). (12) 
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m{ilj,(j)) is called a multiplier and its degree gets subtracted to the degree of the com- 
posed map ^ X 0. A birational map is a rational map with an inverse for this x -product. 
The multiplier for the product of a map and its inverse is the equation of the variety 
where the map is not bijective. When considering iterations of a rational map, we 
always consider this reduced product, which allows for a non-trivial sequence of the 
degrees of the successive iterates. 
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